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c/3 ■ We extend the grid of NLTE models presented in Paper II to explore var i ations 

in abundance pattern in two ways: 1) The adoption of the lAsplund et al\ (|2009[ ) 
(GASSIO) abundances, 2) For stars of metalhcity, [^], of -0.5, the adoption 

^ ■ of a non-solar enhancement of a-elements by +0.3 dex. Moreover, our grid of 

\Q • synthetic spectral energy distributions (SEDs) is interpolated to a finer numerical 

resolution in both Tcs (ATefj = 25 K) and logg (A log (7 = 0.25). We compare 
the values of Tes and logg inferred from fitting LTE and Non-LTE SEDs to 
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ABSTRACT 



O ' observed SEDs throughout the entire visible band, and in an ad hoc "blue" band. 

We compare our spectrophotometrically derived Teg values to a variety of T^s 
calibrations, including more empirical ones, drawn from the literature. For stars 

k>( ■ of solar metallicity, we find that the adoption of the GASSIO abundances lowers 

H ■ the inferred T^s value by 25 - 50 K for late-type giants, and NLTE models 

" " " computed with the GASSIO abundances give Tgfj results that are marginally 

in better agreement with other Te^ calibrations. For stars of [^]=-0.5 there is 
marginal evidence that adoption of a-enhancement further lowers the derived Teg 
value by 50 K. Stellar parameters inferred from fitting NLTE models to SEDs 
are more dependent than LTE models on the wavelength region being fitted, and 
we find that the effect depends on how heavily line blanketed the fitting region 
is, whether the fitting region is to the blue of the Wien peak of the star's SED, 
or both. 

Subject headings: stars: atmospheres, fundamental parameters, late-type 
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Introduction 



We have present e d grid s of LTE (jShort fc Hauschildtl (120101 ). Paper I) and Non-LTE 
(NLTE) ( IShort et all (|2012[ ). Paper II) atmosph eric models and synth etic spectral energy 
distributions (SEDs) computed with PHOENIX fJHauschildt et a/.lll999l ) for GK dwarfs and 
giants (luminosity class V and III) of solar and 1/3 solar metallicity ([^] = 0.0 and —0.5), and 
described a procedure for quantitatively comparing them to mean observed SE Ds carefully 
selected from the large absolute spectrophotometric catalog of lBurnashevI (119851 ) (henceforth, 
B85). Our strongest conclusion from Paper II is that the adoption of NLTE for many opacity 
sources shifts the spectrophotometrically determined T^g scale for giants downward by an 
amount, ATgfj, in the range of about 30 to 90 K all across the mid-G to mid-K spectral class 
range, and across the [^] range from 0.0 to -0.5. This shift brings our spectrophotometrically 
derived Tgfj scale for the solar metallicity G giants into closer agreement with the less model- 
dependent Tgfi scale determined by the IRFM, although our Tes values for these G giants are 
too large in any case. We also found tentative evidence on the basis of two spectral classes in 
the G range that this NLTE downward shift in the T^g scale becomes smaller as luminosity 
class increases from III to V. 



In Papers I and II we divided the observed SED into two spectral bands, designated "blue" 
and "red", with a break-point at 4600 A, which is the wavelength around where the SED 
apparently changes character qualitatively, from being relatively smooth (4600 < A < 7500, 
"red") to being more heavily affected by over-blanketing by spectral lines (3250 < A < 4600, 
"blue"). The quality of the fit of both LTE and NLTE model SEDs to the red band was 
significantly better than that to the blue band. This is to be expected because the blue 
region is more heavily line blanketed, and is expected to be more sensitive to, among other 
things, the distribution of the abundances of individual chemical elements. This region is 
also relatively sensitive to inaccuracies and incompleteness in the input atomic transition 
parameters (oscillator strengths, damping parameters, etc.), and to the treatment of line 
formation and of the atmospheric models {e.g. ID rather than 3D). However, somewhat 
surprisingly, fits to the red and blue bands yielded the same best fit Tcs values for many 
spectral classes. 



In the present investigation we take the step of perturbing the abundance distribution 
of our model grid, by perturbing both the fiducial solar abundance distribution adopted, 
and the choice of a solar or non-solar abundance distribution for the models of [^] = —0.5, 
and studying the effect upon the quality and model parameters of the best fits. We believe 
this is timely and justified as there have been important recent proposed revisions to the 
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solar abundance d i stribu tion , as discussed i n gre ater detail below (for example, compare 
Grevesse fc Sauvall (Il998l ) to lAsplund et al\ ( |2009| )). A goal is to investigate the sensitivity 
of the spectrophotometric Tog scale to the abundance pattern, and to compare it to the 
sensitivity to the thermodynamic treatment (LTE vs. NLTE). This will become increas- 
ingly important as instruments that can produce calibrated spectrophotometry over broad 
wavelength ranges yield data f or larger sarnples of late- type stars {eg. VLT + XSHOOTER 
observations of red supergiants lDavies et al\ ( l2013l )). Furthermore, accurate abundance mea- 
surements in late- type stars, particularly of refractory elements, has become especially impor- 
tant when compa r ing pl anet-hosting stars to the general stellar population (see, for example, 
Adibekyan et al\ ( 120121 ) ). and generally, the accuracy of ab undance measurements de pends 
on the accuracy of the Teg determination (see, for example, iDobrovolskas et al\ (|2012[ ). who 



find from a NLTE analysis that d[|^]/dreff = -0.03 dex/80 K). 



2. Observed /a(A) distributions 



In Paper I we described in detail a procedure for vetting and combining observed SEDs 
selected from the catalog of absolutely calibrated, and uniformly re-calibrated, spectropho- 
tometry (SEDs) of B85, and we only briefly recapitulate here. To state it more starkly 
than we did in Papers I or II, B85 is the largest catalog of absolute stellar spectrophotom- 
etry that we are aware of (~ 1500 stars), and covers a broad enough A range (conserva- 
tively, 3300 — 7500A) to allow evaluation of model SED fits to the heavily line-blanketed 
blue and near-UV s pectr al regions, given a fit to the relatively unblanketed red region. 
Short fc Hauschildtl ( 120091 ) contains a more detailed description of the individual data sources 
included in this compilation. We emphasize again here, that this data set was uniformly re- 
calibrated, by Burnashev (B85), so it should be less affected by random star-to-star errors 
than the original data. Moreover, it contains large enough samples of stars of nominally 
identical spectral type (spectral class and luminosity class), as confirmed by cross-reference 
to other spectroscopic ca talogs, and si r nilar metallicity, as determined by cross-reference to 
the metallicity catalog of ICayrel et al\ ( 1200 ll ). that we can assess the variation, and thus the 
quality, of the SEDs and reject those that are suspect. The final results of our procedure 
are "sample mean" and ±lo" observed SEDs for each spectral type of each metallicity bin 
(A[^] = ±0.1) considered. Paper I contains details of how many SEDs were in each sample, 
and plots showing the typical variation among the observed SEDs, for illustrative spectral 
types. These mean observed SEDs are the "data product" to which we compare our model 
SEDs here. 
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3. Model grid 
3.1. Atmospheric structure and synthetic SED calculations 

3.1.1. Basic grid 

In Paper I we describe an LTE grid of spherical atmospheric models computed with 
PHOENIX for a star of mass IM©, and corresponding high resolution synthetic spectra, that 
spans the spectral class range from about GO to K4 (6125 > Teg > 4000 K, ATgff = 125 K) 
and luminosity class range from V to III (5.0 > logf^ > 1.0, Alogf? = 0.5), and scaled solar 
metallicities, [^], of 0.0 and -0.5, with a mixing length parameter, /, defined by the Bohm- 
Vitense treatment of convection of one pressure scale height. In the current investigation, 
we are concerned with the sub-set of the grid corresponding to giants (3.0 > \ogg > 1.0) 
(the only metal-poor ([^] ^ —0.5) GK stars in the B85 catalog are giants and supergiants, 
and the modeling of supergiants requires special additional considerations that are beyond 
the scope of this investigation.) The choice of [^] values was guided by the availability 
of objects in the B85 catalog. Paper II describes the NLTE version of this grid, in which 
the lowest two ionization stages of 24 chemical elements up to Ni, and thousands of atomic 
transitions {b — b and b — f) are treated in self-consistent NLTE statistical equilibrium (SE), 
accounting simultaneously for NLTE radiative transfer in those radiative atomic transitions 
that are accounted for in the SE rate equations. 



3.1.2. Interpolated SED grid 

In Paper I, we described a finer grid of SEDs, with AT^s = 62.5 K, derived from 
this basic grid by linear interpolation in the log/A(Teff) relation. Here we have derived an 
even finer grid, with ATeg = 25 K and A log (7 = 0.25, by quadratic interpolation in the 
log/A(Teff) and linear interpolation in the /A(logf7) relations. The log/A(Teff) interpolation 
was performed using two different methods to interpolate values in the interval [yi,yi+i]: i) 
Simple quadratic interpolation by fitting a parabola to the three points i/i-i, i/i, Vi+i, and 
ii) Least-squares quadratic interpolation by fitting a parabola to the four points t/j-i, i/i, 
Vi+i, yi+2- The two methods generally gave interpolated log /a values that were within 1% 
of each other, and were both found to yield an interpolated log /a (Teg) function that was 
smooth and well-behaved as determined by exhaustive visual inspection. Moreover, SEDs 
interpolated with both methods yielded the same best fit parameters to within the numerical 
precision of our grid when the x^ deviation from observed SEDs was minimized (see Section 
H]). The results presented here are based on the least-squares quadratic interpolation. 
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The choice of AT pff = 25 K is consistent with the recent ground-breaking investigation of 
Huber. et all ( 120121 ). who found mean deviations of 22 ± 32 K among values of Tefj derived 
with interferometric-asteroseismic, spectroscopic, and photometric techniques for stars of 
4600 < Teff < 6200 K, including red giants. 



Synthetic SEDs and microturbulence Synthetic SEDs suitable for comparison to the 
observed SEDs are prepared by convolution with a Gaussian kernel with a FWHM value 
of 75 A, which was found to be the value, to within ±5 A, that produced synthetic SEDs 
most closely matching the appearance of the observed SEDs. The nominal sampling of the 
B85 catalog is 25 A, but the instrumental profile is apparently unknown. In Paper I and II 
we presented a grid in which the microturbulent broadening parameter, ^t, increases from 
2.0 to 4.0 km s~^ as log*? decreases from 3.0 to 1.0. To remove variation of a physically 
ill-motivated variable from our grid, we have re-computed the atmospheric models and the 
synthetic spectra in this logg range with ^t = 2.0 km s~^. Our tests show that for SEDs 
broadened for the observed spectral resolution element, AA, of 75 A, the difference between 
a SED computed with C,t = 2.0 and 4.0 km s~^ is completely negligible. 

Finally, both the observed and smoothed synthetic spectra are put onto the same flux scale 
by whole-area normalization following the procedure described in Paper II. For illustrative 
purposes. Fig. [3] shows a plot of the interpolated log/A(log A,Teff) surface, broadened to 
match the spectral resolution of the data and normalized, as described above, for logg = 2.0 
for NLTE models with the GS98 abundances. The coarser grid from which this finer grid was 
interpolated is also shown for reference. We show the broadened log f\ surface because the 
unbroadened surface is highly affected by spectral line blanketing, and is visually confusing. 



3.2. Abundances 



Solar distributio ns The grid presented in Papers I and II was computed using the 
solar abundances of iGrevesse fc Sauvall (Il998[ ) (GS98, henceforth), and this distribution was 
simply scaled for the [^] = —0.5 models. The GS98 distribution is based on a compilation 
and review of abundances derived largely from ID LTE semi-empirical solar photosphere 
and spectrum modeling, checked against meteoritic abundances for relevant elements. These 
abundances, along with other similar distributions published by the same authors over the 
years, have be en widely used . We have re-co mputed our [w ] = 0. grid with the recent 



abundances of lAsplund et all ( l2009l ) (see also iGrevesse et all (120101 )) (hereafter GASSIO), 
which are based on recent 3D hydrodynamic solar photosphere and spectrum modeling, and 
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which account for NLTE for some species, particularly C and O. The GASSIO distribution 
has lower abundances, by 0.1 dex or more, than that of GS98 for the important elements 
C, N, O, Ne, Na, S, and K. In particular, Na is among the most important e~ donors in the 
solar atmosphere. Given the scale of the computational task of re-computing ~ 500 models 
and spectra, we have restricted ourselves to investigating the effects of only these two solar 
distributions, and take them as being representative of the recent variation in the reputable 
measured solar abundance distributions. 

Fig. [T] shows a comparison, with residuals, of the area-normalized, smoothed NLTE SEDs 
computed with the GS98 and GASSIO abundances for solar abundance models of log (7 
equal to 2.0 and T^s equal to 4000, 4625, and 5250 K {ie. spanning almost the entire T^s 
range at a middle value of log^?). (Note that the residuals vary with large amplitude in A 
regions where the flux derivative -^ of the SED is large.) The main effect of perturbing 
the abundance distribution is that the GASSIO abundances lead to a larger blue and near- 
UV band flux as compared that of the GS98 abundances, while the two distributions yield 
similar red band fluxes. This effect becomes more pronounced as Teff decreases. This is 
to be expected given that the GASSIO distribution generally has lower abundances than 
the GS98 distribution for many relatively abundance "light metals" that contribute to the 
line blanketing, including Ti, V (both A[f ] = -0.07), Cr, and Ni (both A[f ] = -0.03), 
thereby reducing the line extinction in the most heavily blanketed spectral region (the blue 
and near-UV), thus allowing more flux to escape there. This is essentially a reduction in 
the classical line-blanketing/backwarming effect in the GASSIO models as compared to the 
GS98 models. Therefore, we expect that Teg values inferred from fitting GASSIO models will 
be lower than those of GS98 models to compensate for GASSIO models having a "warmer" 
SED at a given model Teg value. 



Metal poor models For our grid of [^] = —0.5 models, we investigate the effects of four 
different abundance distributions: simply scaled solar distributions adopting the GS98 and 
GASSIO abundances, and a-enhanced distributions in which the scaled distributions of GS98 
and GASSIO are altered to a non-solar distribution by enhancing the abundance of eight a- 
process elements (O, Ne, Mg, Si, S, Ar, Ca, and Ti) by 0.3 dex. We designate the latter 
two distributions GS98-a and GASSlO-a. Chemical abundance analyses of kinematically 
selected halo and thick disk stars in the solar neighborhood, including giants, have found 
that stars of [^1 < —0.5, can h a ye g-element abu n dance s enhanced by as mu c h as ~ 0.4 



dex (s ee lAdibekvan et al\ J20ia). lAdibekvan et al\ J2012h . iRuchti et a/l feoilj ). iFuhrmann 



Furthermore, 



( 2OIII ) and lNissen fc Schuster] (120101) . for examples based on large surveys and LTE analyses). 



Peterson et al\ ( 119931 ) performed a detailed LTE spectrum synthesis analysis 



-7- 



of the high resolution spectrum of the bright standard star Arcturus ( a Boo), of Jf|-] ~ — 0.5 



and found g-element abundanc es enhanced by ~ 0.3 dex. Recently, IShi et all ( 120121 ) and 



Takeda fc Takada-Hidai I (120121 ) have carried out NLTE analyses of the abundance of the a- 
elements Si and S, respectively, and confirmed similar enrichment in moderately metal poor 
stars. We have chosen an a-enhancement at the high end of the range that has been found 
for moderately metal-poor solar neighborhood stars (+0.3 dex), as exemplified by Arcturus 
and applied it to all elements of even atomic number from O to Ti, to deliberately maximize 
any effect on derived stellar parameters caused by a-enhancement. 



Fig. |2] shows a comparison, with residuals, of NLTE SEDs computed with the GS98 and 
GASSIO abundances for metal-poor models at a logg value of 2.0, and Tes values of 4000 
and 5000 K (similarly to Fig. [1]). The blue/near- UV brightening caused by adoption of 
the generally lower GASSIO abundances is again evident, despite the reduced significance 
of line extinction in moderately metal-poor models. We expect the same reduction in the 
inferred best fit T^s value as for the solar abundance case, but with reduced magnitude. For 
the GS98 models, we also show the SED with and without a-enhancement. The effect of 
a-enhancement is to decrease the blue/near- UV band flux with respect to the red band flux, 
yielding "cooler" SEDs for a given model T^e value. We expect a-enhancement to have the 
same effect as changing from the lower GASSIO abundances to the higher GS98 abundances 
in that it should lead to "cooler" SEDs for a given model Teff and therefore, a larger best-fit 
Teff value. 



4. Goodness of fit statistics 

We minimize x^ fo^' the difference between the observed SED and the trial synthetic 
SEDs (the "hypotheses"), relative to the observed SEDs (ze. (/A,obs — /A,Mod)//A,obs) in each 
wavelength sampling element, AA. This is equivalent to adopting unity {f\,ohs/f\,ohs) as the 
"observed" quantity, and (/A,Mod//A,obs) as the trial quantity; ie. 



x' 



1 ^ 

^^ = ir7 2^((/A,Obs — /A,Mod)//A,Obs) /(/A,Mod//A,Obs) (1) 



A 



or, equivalently. 



X ^^ — T7 2^(1-0 ~ /A,Mod//A,Obs))/(/A,Mod//A,Obs) (2) 



A 



where N is the number of wavelength elements minus 1, which we take to be the number 
of degrees of freedom. The sampling interval, AA is 15 A, which was chosen to slightly over- 
sample SEDs broadened to match a spectral resolution element of 75 A (see Section 13.1. 2p . 
For each trial model we compute three values of x^AA~^, one each for our "blue" band (3250 
- 4600 A), our "red" band (4600 - 7500 A), and the entire range, to study the quality of the 
fit to the over-blanketed blue region with respect to the red region. For fits to the overall 
spectrum, the minimum value of x^^-^~^ ranged from 0.002 to 0.005, with the exception 
of the coolest, most line blanketed sample, K3-4 III, for which the x^^-^"^ minimum was 
Ri 0.01. These values are also characteristic of the fit to the "blue" band alone, indicating 
that the quality of the fit to the more heavily line-blanketed blue band dominates the quality 
of the global fit. 

Paper II contains a detailed comparison of the stellar parameters derived from fitting 
NLTE model SEDs as compared to LTE SEDs, and we do not repeat that discussion here, 
except to note that the main result was that NLTE T^s values derived from fitting the 
overall visible band were found to be systematically lower than LTE values by one numerical 
resolution element, AT^s, in the coarser SED grid of Paper I (nominally 62.5 K). With the 
finer numerical resolution of our current SED grid, we find that this Tes difference for the 
GS98 models generally varies from 50 to 75 K (in the new grid, this amounts to two to three 
numerical resolution elements, ATcfj). This is a systematic offset between the NLTE and 
LTE SED-derived T^s scales that we are now numerically resolving, whereas in Paper II it 
was only marginally resolved. 

Figs, m and [5] show the logarithm of the x^AA^^ value as a 2D surface plotted versus 
model Tes and log g in the vicinity of the x^ minimum for the G5 III and K3-4 III samples 
of [^] = 0.0, respectively, for NLTE models with the GS98 abundances. We have chosen a 
logarithmic scale to enhance the relatively weak variation of x^ with log g with respect to its 
variation with T^fi. Visual inspection of such plots for all our spectral class samples assures 
us that the x^ minimum is global and well defined as a function of Tgfj. The minimum is less 
well-defined as a function of logg, and in every case we have checked for local degeneracy in 
Teff /logs' around the x^ minimum (see the discussion of the Kl III/[^] = 0.0 sample in this 
section) . 

Table [1] shows the best fit Teg and logs' values for the models of [^] = 0.0 for the GS98 
and GASSIO abundance distributions in both LTE and NLTE. Figs. [6] and [7| show the 
variation of x^AA~^ with NLTE model T^s and logs, respectively, for the stars of solar [^] 
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value for the models of GS98 and GASSIO abundance distributions, and Fig. [8] shows the 
best fit NLTE T^s value as a function oi B — V color for these stars. Generally, the NLTE 
GASSIO abundances pervasively yield best fit T^s values that are 25 to 50 K (one or two 
ATefj resolution elements) cooler than the older GS98 abundances, while yielding the same 
best fit log g values, a difference that we are marginally resolving. We conclude that there is 
evidence that adoption of the GASSIO abundances leads to best fit NLTE Tgg values for the 
blue band and total SEDs that are 25 to 50 K cooler as compared to those with the GS98 
abundances. This is consistent with what we expect from the comparison of GASSIO and 
GS98 SEDs discussed in Section 13.21 Generally, the two abundance distributions yield the 
same NLTE best fit log*? value to within the A log (7 grid resolution of 0.25. 

The one exception is the Kl III sample, for which the GASSIO models yield a Tcs value 
that is 50 K warmer, while yielding a logg value that is 0.5 smaller (two Alogf? resolution 
elements). We note that the best fit logg value of the GS98 models (2.50) is higher than 
those of the other spectral classes, for either abundance distribution, and that the value of 
X^ for the Tes/logg pair 4600/2.0 is identical to three digits after the decimal with that for 
the 4550/2.5 pair. If the logg value is constrained to be 2.0 (GASSIO best fit model value) 
then the NLTE GS98 models yield a Tefj value of 4600 K, in closer agreement with the GS98 
vs GASSIO behavior of the other spectral classes. 

From Fig. [HI which also includes best fit LTE Tcs values, we note that the effect of 
NLTE on derived Tes values compared to LTE is about the same for both the GASSIO and 
GS98 abundance distributions, namely a reduction of 50 to 75 K. This is to be expected 
because the GS98 and GASSIO distributions have identical values for the Fe abundance, and 
a brightening of the SED in the blue band with respect to the red band c aused by NLTE 



over-ionization of F e I is the dominant NLTE effect in late-type stars (see iRutteru ( 1l986l ) 



Short fc Hauschildtl tom ). 



[^] = —0.5 Table [2] shows the best fit T^q and logg values for the models of [^] = —0.5 for 
all four relevant abundance distributions in both LTE and NLTE. Fig. [9] shows the variation 
of x^AA~^ with NLTE model Tgfj for the stars of [^] equal to -0.5 for models of scaled GS98 
and GASSIO distributions and the GS98-a and GASSlO-a distributions, and Fig. [TOl shows 
the best fit T^g values as a function oi B — V for both LTE and NLTE modeling. For the 
G8 sample, adoption of the GASSIO abundances decreases the derived T^s value by 25 K, 
similarly to its effect on the [^] = 0.0 models, whereas for the K1.5 III sample, the GASSIO 
abundances increase the Tgs value by 25 K. 
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ck enhancement More clearly, for the G8 sample, the effect of a enhancement is con- 
sistently to further lower the inferred Teg value, by 50 K for both the GS98 and GASSIO 
abundances. By contrast, the effect of a enhancement on the K1.5 III sample is less clear; 
for the GS98 models it increases the derived Tgff value by 25 K, whereas its affect of the 
GASSIO models is to lower T^s value by 25 K. Given that one numerical resolution element, 
ATeff, is 25 K, we conclude that for stars of [^] = —0.5, we are not resolving a clear trend 
for either the effect of choice of solar abundance distribution, or of a enhancement, with the 
possible exception of the effect of the latter on the derived T^q values of G8 III stars. 

The choice of solar abundance distribution has a negligible effect on the derived \ogg 
values. However, there is evidence that a-enhanced models lead to derived \ogg values that 
are 0.25 to 0.5 (one to two A log (7 resolution elements) larger than those of scaled solar 
models, regardless of choice of solar abundance distribution. 

As with the models of [^] = 0.0, the NLTE reduction of the Teg scale is independent of 
choice of solar abundance distribution, although we note that the effect is more pronounced 
for the G8 III sample (75 K for three of the four abundance distributions), than it is for the 
K1.5 III sample (25 K for three of the four abundance distributions, and never more than 
50 K). 

Blue band Fitting the "blue" band with LTE GS98 and GASSIO models results in best 
fit Tcff and \ogg values that are generally the same as those from fitting the overall visible 
band, with occasional deviations by one numerical resolution element (25 K and 0.25, re- 
spectively) higher or lower. The "blue" band and the overall visible band yield the same 
stellar parameters to within our ability to resolve differences. This is not surprising because 
the heavily line blanketed "blue" band is most sensitive to Tgff, and it thus dominates the 
quality of the global fit. For the NLTE models, GS98 and GASSIO models fit to the "blue" 
band yield best fit Teg values that are larger by 25 to 75 K (one to four resolution elements) 
than those fit to the overall visible band for two of the six spectral class samples (G5 and 
Kl). The larger "blue" band T^s value tends to be accompanied by a best fit \ogg value that 
is lower by 0.25 to 0.5. This result applies to samples of both metallicities. We conclude 
that stellar parameters inferred from fitting NLTE models to SEDs may be more dependent 
than LTE models on the wavelength region being fitted, and find that the effect depends on 
how heavily line blanketed the fitting region is, whether the fitting region is to the blue of 
the Wien peak of the star's SED, or both. 
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5. Comparison to other Teg calibrations 



Figs. [8] and [TT] show the comparison of our derived T^s values with those of other 
investigators for the solar metallicity stars. Papers I and II contains a review of the recent Tgs 
calibrations in the literature, based on a var iety of complementary techn iques, to which we 
compare our results. To recapitulate briefly, iRamirez &: MelendezI (120051 ) present a Teg (5 — 
V) calibrati o n for a wide range of stars determined with the Infrared Flux Method (IRFM), 
Wang et al\ (120111 ) present Tes values for G giants determined fr o m bot h B — V and Fe I/II 
ionization equilibrium along with B — V values, iBaines et all (J2010l ) present T^s values 
for select K giants d etermined from optical interferometry along with MK spectral classes, 
Takeda. et al\ (120081 ) present T^s values fo r G giants determined from modeling Fe I and II 
lines for stars of given spectral class, and JMishenina et al\ (120061 ) present Tes values for G 
giants determined from line depth ratios for stars of given spectral class. We note that Paper 
II contains a thorough investigation of how well our NLTE T^s values compare to the results 
of other studies with respect to our LTE T^s values. Although we include the LTE results 
again here in Figs. [S] and [TT] for interest, we do not repeat the discussion of NLTE versus 
LTE. 



To facilitate the comparison to the results of iRamirez &: MelendezI (120051 ) and lWang et al 
(I2OIII ) in Fig. [HI we determine the me an observed B — V color of our spectral class samples 
using the catalog of iMermilliodI (1199 ll ) following the procedure described in Paper I. In Fig. 
[HI among the NLTE models, for every class except K2 III, NLTE models with the GASSIO 
abundance distributi on yield Tpff values that are c loses t to, or are among the closest to, the 
values determined by [Ramirez fc MelendezI (|2005| ) and IWang et al\ (1201 ih . 



From Fig. [TT], among the NLTE models, for two of our three spectral classes that overlap 
with the interferometric results of iBaines et al\ (120101 ). at the cool end of our spectral class 
range, the GS98 abundances lead to NLTE T^^ values that are in better agreement with that 
study. For the G5 stars, at the other end of our ra nge, the NLTE GASS IO ab undances lead to 



Tpff va lues that are closer to the G giant results of iTakeda. et al\ (120081 ) and iMishenina et al. 



(l2006l ). although any conclusion is weakened because all our T^s values are generally too 
high compared to other results at the warm end of our Teg range. 



From Fig. [TD] we can see that for the K1.5 sample, which consists entirely of multiple 
independent observations of Arcturus (a Boo, see Paper I), among the NLTE models, the 
G ASSIO and GS98-Q; di s tribut ions give a Teg value that is closest to the IRFM calibration 
of [Ramirez fc Melendea (|2005[ ). That the GASSlO-a distribution yields a more discrepant 
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Tpff value is no t e-wor thy in that a Boo has been found to be a-enhanced (see, for example, 



Peterson et all Jl993h ) 



6. Conclusions 

As expected, for solar metallicity stars, we find that the effect of NLTE on spectropho- 
tometrically derived Teg values compared to LTE is about the same for both the GASSIO 
and GS98 abundance distributions, namely a reduction of 50 to 75 K. As with the models 
of [^] = 0.0, for the metal-poor stars the NLTE reduction of the Tgfj scale is independent of 
choice of solar abundance distribution, although we note that the effect is more pronounced 
for the G8 III sample (75 K for three of the four abundance distributions), than it is for the 
K1.5 III sample (25 K for three of the four abundance distributions, and never more than 
50 K). 

For solar metallicity stars, we conclude that there is evidence that adoption of the GASSIO 
abundances lowers the inferred NLTE spectrophotometric T^g scale by 25 to 50 K as com- 
pared to the GS98 abundances. Generally, the two abundance distributions yield the same 
NLTE best fit logg value to within the Alogf^ grid resolution of 0.25. For solar metal- 
licity stars, NLTE models computed with the GASSIO abundances give, on balance, T^s 
results that are marginally in better agreement with other Tcfj calibrations in the literature 
determined by a variety of means. 

We conclude that for stars of [^] = —0.5, we are not resolving a clear trend for the effect 
of choice of solar abundance distribution on the Tes value inferred from the SED. There is 
marginal evidence, based entirely on the G8 III sample, that the adoption of a-enhancement 
leads to a reduction of 50 K in the NLTE SED-fitted T^g value. The choice of solar abundance 
distribution has a negligible effect on the derived logf^ values. However, there is evidence 
that a-enhanced models lead to derived logg values that are 0.25 to 0.5 (one to two Alogf? 
resolution elements) larger than those of scaled solar models, regardless of choice of solar 
abundance distribution. With only two spectral classes, it is difficult to draw conclusions 
about which abundance distribution better fits other Teg calibrations. 

We conclude that stellar parameters inferred from fitting NLTE models to SEDs are more 
dependent than LTE models on the wavelength region being fitted, and find that the effect 
depends on how heavily line blanketed the fitting region is, whether the fitting region is to 
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the blue of the Wien peak of the star's SED, or both. 



We will make both the library of observed SEDS and the NLTE (and corresponding LTE) 
grid of model SEDS with all abundance distributions available to the community by ftp 



(http://www.ap.smu.ca/ ishort/PHOENIX). 



CIS is grateful for NSERC Discovery Program grant 264515-07. The calculations were 
performed with the facilities of the Atlantic Computational Excellence Network (ACEnet). 
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Fig. 1. — Comparison of normalized, smoothed (see text) SEDs for the 3200 to 4500 A region 
computed with different abundance distributions for models of [^] = 0.0 and logg = 2.0. 
Thin lines: GS98 models; thick lines: GASSIO models; dark lines: T^s = 4000 K; interme- 
diate lines: T^s = 4625 K, lightest lines: T^s = 5250 K. Upper panel: Normalized SEDs 
(see text for description of normalization); lower panel: residuals (relative differences) with 
respect to the GS98 models. Note that the residuals show large variation in regions where 
the flux derivative, -^, is large. 
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Fig. 2. — Same as Fig. [H except for models of [^] = —0.5, at Teg = 4000 K (dark lines) 
and 5000 K (light lines). For the GS98 models, we also include SEDs computed with a- 
enhancement (dashed lines). In the lower panel all residuals are with respect to the GS98 
models without a-enhancement. 



-17- 



0-6 



0.0 



^5 c 



Q 



^1.0 



^1-S 




^■'^'";\^'■'^■■'\■■^■■'■'*■«^jv■^■>^^■-^■'/•^•V■^''■V■■V•^■'>^ 
^"^■^■^;'■'•■■'■:^v;;:..^^^;;^^:;^^\■-;;^;-■<■^■;^v^■^■^>^,;■.:•.',-.>^:i■■;,■. 

^;,;^;,^:^:.;;,%-;;;;3^^.^;■::.;■::;,:;:,■V,^;yO 

•.'.■'■■ v. •■'.'.d....,wii '.■■'■. !■■.■-.' '.>.'■*:;. ■'■.■.'■". '.'■.'■■* " 



,•.■..■•••■■■■"'■«•:■.,■ 
:^,:::;!::::m. 



^^-■■■■'•'/-■.•/■■.^■•^■j..' .■■■■:■.■:■■:■■:■■■ 
■■:■ ■ ■■::>;/'>^N%N'::>: 



<=>^' 



.oo 




^O" 



>o^^^7^^ 



so*^ 



Fig. 3. — The normalized, smoothed (see text) log /A(log A) surface as a function of log A and 
model Teff for the fine grid of solar abundance NLTE models with the abundances of GS98 
and logg = 2.0 interpolated to ATcs = 25 K, for illustrative purposes. The coarser basic 
grid among which the finer grid was interpolated (ATcfj = 125 K) is also shown (thicker 
lines) . 
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Fig. 4. — The logx^AA^^ surface as a function of trial Tes and logg values in the vicinity 
of the x^ minimum for the G5 III/[^] = 0.0 sample and NLTE models of GS98 abundances. 
Loci of constant T^s and log^f are shown with solid and dotted lines, respectively, and the 
model Teff /log g pair of minimal x^ value is marked with a '+' symbol. The loci of of constant 
Tgff and log g corresponding to the best fit log g and Tcfj values, respectively, are plotted with 
thicker lines. We plot logx^ to enhance the weaker variation of x^ with model logg as 
compared to its variation with T^s. 
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Fig. 5.— The same as Fig. Shut for the K3-4 III/[f ] = 0.0 sample. 
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Fig. 6. — The x^AA~^ value as a function of trial Tcs value for stars of [^] = 0.0. We 
show x^ for NLTE models with the GS98 and GASSIO abundances (solid and dashed lines, 
respectively). The results for LTE models with the GS98 abundances are also shown for 
comparison (dotted line). 
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Fig. 7. — Same as Fig. |6l but showing "X^lSX ^ as a function of log (?. 



-22- 



0) 



5500 


GS98 ' 
"GA10 


2.25 2.00 
2.25 2.00 


1 ' 

1.50 
1.50 


1 1 1 1 

2.502.00 
2.002.00 


1 1 1 1 

1.75 
1.75 




^ 








- 


5000 


* 


v □ 






- 




^ 


^ 






- 




■ 






. O 


■ 




■ 






^s° 


■ 


4bUU 


■ 








- 




■ 


G5III 


KOIII 


K2III 


No 


4000 


1 


G8III 


1 


K1III 

1 


1 



0.8 



1.0 1.2 

B-V 



1.4 



Fig. 8. — Best fit Teg values for stars of [^] = 0.0 from this investigation for NLTE (light 
gray) and LTE (dark gray) models with the GS98 (squares) and GASSIO (diamonds) solar 
abundances. For comparison we sh ow the results of othe r inves tigations of the Tefj (i? — V) 
calibration: the IRF M ca libration of lRamirez fc MelendezI ( 120051 ) (solid black line) and values 



from IWang et al\ (l201ll ) derived from B — V photometry (triangles) and from Fe I/Fe II 



ionization equilibrium (asterisks). Our spectral classes have been assigned corresponding 
B — V values following the procedure described in Paper I to enable comparison with the other 
studies. The annotation shows the best fit \ogg values derived from each solar abundance 
distribution ("GAIO" indicated GASSIO abundances). 
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Fig. 9. — Same as Fig. [6] but for stars of [^] = —0.5. We have included results for NLTE 
models with the a-enhanced GASSlO-a abundances (dot-dashed lines). For clarity, we have 
omitted the NLTE models with GS98-a abundances, however, the complete results in Table 
[1] show that the effect of a-enhancement is similar for the GS98 and GASSIO models. As 
with Figs. [6] and [71 we include the LTE models with GS98 abundances for comparison. 
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Fig. 10. — Same as Fig. |H1 but for stars of [^] = —0.5. Results derived from the GS98-a 
and GAS SlO-g models are denoted w ith plus signs ('+') and 'X' symbols, respectively. The 
results of lRamirez &: MelendezI ( 120051 ) and I Wang et al\ ( 120111 ) are those for the corresponding 
metallicity. 
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Fig. 11. — Same as Fig. [HJ but with a comparison to the res ults of other invest igations 
of the Tpff (spect ral class) cahb r ation : Interferometric results of iBaines et ali (120101 ) (trian- 
gles), results of JTakeda. et al\ (120081 ) from Fe l/Fe II ionization equilibrium (X's), and of 



Mishenina et al. 



( 20061) frorn line depth ratios (asterisks). For the results of IXakeda. et al. 



(120081 ) and iMishenina et al\ (120061 ) we also show the mean T^^ value for each spectral class 
(large, dark X's and asterisks, respectively). 
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Table 1. Best fit parameters {Tes/logg) for total and "blue"-band SEDs: Models of 

[M] = 0.0 







LTE 






NLTE 






GS98 


GASSIO 


GS98 


GASSIO 


SED: 


Total 


Blue 


Total 


Blue 


Total 


Blue 


Total 


Blue 


G5 


5350/2.25 


5375/2.25 


5325/2.25 


5325/2.25 


5300/2.25 


5325/2.00 


5250/2.25 


5275/2.25 


G8 


5150/2.00 


5150/2.00 


5100/2.25 


5075/2.25 


5075/2.00 


5075/2.00 


5025/2.00 


5025/2.00 


KO 


4900/1.50 


4900/1.50 


4850/1.75 


4850/1.75 


4825/1.50 


4825/1.50 


4800/1.50 


4800/1.50 


Kl 


4625/2.50 


4575/3.00 


4675/2.00 


4700/1.75 


4550/2.50 


4625/2.00 


4600/2.00 


4675/1.50 


K2 


4575/1.75 


4575/1.75 


4525/2.00 


4525/2.00 


4475/2.00 


4475/2.00 


4450/2.00 


4450/2.00 


K3-4 


4250/1.75 


4250/1.75 


4175/2.00 


4175/2.00 


4200/1.75 


4200/1.75 


4150/1.75 


4125/2.00 





















Table 2. Same as Tabled! but for total SED fits only: Models of [-^] = —0.05 



LTE 



GS98 
Solar 



a+ 



GASSIO 
Solar 



a+ 



GS98 
Solar 



NLTE 



a+ 



GASSIO 
Solar 



a+ 



G8 4750/2.00 4675/2.50 4725/2.00 4650/2.50 4675/2.00 4625/2.50 4650/2.00 4600/2.25 
K1.5 4275/1.50 4325/1.25 4350/1.00 4300/1.25 4275/1.25 4300/1.25 4300/1.00 4275/1.25 



